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Abstract 



We consider a gauge/gravity dual model of thermalization which consists of a col- 
■ lapsing thin matter shell in asymptotically Anti-de Sitter space. A central aspect of 

qq ■ our model is to consider a shell moving at finite velocity as determined by its equa- 



tion of motion, rather than a quasi-static approximation as considered previously in 
the literature. By applying a divergence matching method, we obtain the evolution of 
singularities in the retarded unequal time correlator G R (t,t'), which probes different 
stages of the thermalization. We find that the number of singularities decreases from 
a finite number to zero as the gauge theory thermalizes. This may be interpreted as a 
sign of decoherence. Moreover, in a second part of the paper, we show explicitly that 



the thermal correlator is characterized by the existence of singularities in the complex 
time plane. By studying a quasi-static state, we show the singularities at real times 
originate from contributions of normal modes. We also investigate the possibility of 
obtaining complex singularities from contributions of quasi- normal modes. 
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1 Introduction and Summary 



Recent experimental results have shown that collective phenomena observed in heavy ion 
collisions at the relativistic heavy ion collider (RHIC) are well-described by hydrodynamics. 
The use of hydrodynamics assumes a short thermalization time, r ~ 0.5/m, which is the 
time scale for the matter produced in the collisions to reach local equilibrium. The new 
experiments using heavy ion collisions at the Large Hadron Collider (LHC) are expected to 
provide further constraints on the thermalization time. The short thermalization time is 
believed to be due to the strongly coupled nature of the matter produced in the collisions. A 
theoretical understanding of the thermalization mechanism requires knowledge of Quantum 
Chromodynamics (QCD) dynamics at strong coupling and far from equilibrium, which is in- 
accessible to lattice simulations and to perturbative field theory techniques. Gauge/gravity 
duality offers a useful tool to study the dynamics of strongly coupled gauge theory. Within 
the framework of gauge/gravity duality, the gravity dual of N = 4 Super Yang- Mills theory 
at finite temperature is given by the AdS-Schwarzschild black hole geometry. It is there- 
fore natural to conjecture that the thermalization process is dual to black hole formation 
via gravitational collapse in Anti-de Sitter space. Solving the necessary Einstein equations 
generally requires numerical computations. Work along this line includes [H EJ (3J IU [U EJ [7] . 
While these works provide valuable information on the evolution of the one-point function 
in the gauge theory undergoing thermalization, perhaps equally important is the evolution 
of the two-point function, which encodes information on the correlation and spectrum of a 
given operator. The study of the two-point function amounts to a study of the behavior 
of a bulk field in a gravitational collapse background. Recently, there have been intensive 
efforts in studying the evolution of the two-point function within gauge/gravity duality. 
An incomplete list in the context of thermalization in heavy ion collisions can be found in 
PElIinillllliaiiaillliailSIiailElIiaEO]. In particular, equal-time quantities have been 
extensively studied in [15], see also [20], which showed interesting patterns of spatial corre- 
lation. On the other hand, unequal-time correlators, in which the operators are inserted at 
different times, allow to probe the causal structure and temporal correlation of the gauge 
theory and are therefore complementary to equal-time quantities. In [16j . an initial value 
problem was formulated for the unequal-time correlator and was applied to the background 
of an evolving black hole [21] . 

In this paper, we study the gravitational collapse of a massive shell as in [11] . We 
include the finite velocity motion of the shell in our analysis, thus going beyond the quasi- 
static approximation used previously. The study of a massive shell has the advantage 



2 



that the corresponding field theory correlators have manifest singularities. This feature 
is absent in other collapsing shell models, for instance for the Vaidya metric, where the 
shell is massless. For studying the singularities in the gravitational collapse model, we use 
techniques which we developed in our previous papers, where we considered the toy model of 
a mirror moving in Anti-de Sitter space. In [13], we have calculated the spatially-integrated 
unequal time correlator by studying a bulk scalar in a moving "mirror" background, that is 
the scalar satisfies the Dirichlet boundary condition on a prescribed surface. We found that 
in the WKB limit, the singularities of the resulting correlator showed a pattern consistent 
with a geometric optics picture: The singularities occur at those times a light ray bouncing 
between the AdS boundary and shell is reflected at the boundary. For a static mirror, this 
was already realized in [22]. The relation between singularity locations and geometric optics 
can be viewed as a realization of the bulk-cone singularity conjecture formulated in |10] , 
In a follow-up paper p3], the present authors and Hoyos developed a powerful divergence 
matching method, which allows for the determination of the precise form of the singularities 
without solving the equation of motion for the scalar in the bulk. This has paved the road 
for studying of the structure of the singularities in a gravitational collapse model. 

In this paper, we generalize the divergence matching method in order to apply it to 
the case of gravitational collapse. We focus on the singularities of the retarded correlator 
Gn(t, t') = 9(t — t')([0(t), 0(t')]). For varying t' , we analyze the singularities of the correla- 
tor Gji(t, t') and find that singularities occur at t = i n , where i n (ii < I2 < ■ ■ ■ ) is the time 
for a light ray originally leaving the boundary at t' to return to the boundary after the n-th 
bouncing off the shell. The main result we find using the divergence matching method is an 
evolution of the singularities in the process of thermalization. We will see that unlike in the 
problem of a moving "mirror" background, where there are infinite number of singularities 
in the correlator, the case of gravitational collapse only contains a finite number of singular- 
ities. In particular, the n-th singularity of i n moves monotonously to +00 as t' approaches 
a critical value T n (Ti > T2 > • • • ) from below. For t' > T±, the last singularity ii escapes 
from detection and the correlator appears thermal as far as the singularities are concerned. 
Since the singularities in the unequal time correlator measure the strongest correlation in 
time, the disappearance of singularities has the interpretation of temporal decoherence. 

A second aspect we study in this paper is the correlator for a field in thermal equilib- 
rium, which is characterized by the appearance of singularities in the complex t plane, which 
is closely related to the structure of the quasi-normal modes (QNM). The geometric optics 
picture has been used in [22] to extract the asymptotic QNM in the AdS-Schwarzschild 
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background, where the asymptotic QNM are obtained as the reciprocal of the complex time 
period of a light ray bouncing in AdS-Schwarzschild background. We confirm this picture by 
explicitly evaluating the retarded correlator and identifying the singularities in the complex 
t plane. To understand the appearance of the singularities in the complex t plane, we further 
investigate the evolution of the QNM for the quasi-static states dual to a shell levitating 
at rest at different positions above the horizon. We study the evolution of the QNM as the 
shell is lowered to the horizon. We find that among the QNM the normal modes allows us 
to reproduce the singularities for real t obtained by the divergence matching method, while 
the complex QNM does not lead to the expected singularities in the complex t plane. 

The paper is structured as follows: In Section 2, we develop a differential form of the 
divergence matching method, which is more suitable for generalization to the gravitational 
collapse model. In Section 3, we review the gravitational collapse model used in and 
study the behavior of a bouncing light ray in the collapsing background. Then we generalize 
the divergence matching method to the case of gravitational collapse model and check that 
it passes a non-trivial test. The application of the method gives rise to the singularity 
evolution described above. Section 4 is devoted to the calculation of the unequal-time 
correlator in thermal equilibrium, which shows the appearance of the singularity in the 
complex time plane. In Section 5, we revisit the results of Section 3 and give a dual picture 
of the evolution of the QNM with an explicit example of a sequence of quasi-static states. We 
discuss open questions and future directions in Section 6. Some notes on the computation 
of black hole QNM are collected in the appendix. 

2 Differential form of divergence matching for Dirichlet prob- 



Let us first reformulate the divergence matching method developed in |14| for the Dirichlet 
problem in a differential form, which is easily generalizable to more complicated models. 
We start with the Dirichlet problem 



Here z = f(t) is the Dirichlet surface, which can be viewed as a mirror. Without the mirror, 



lem 




(i) 



G R (t,z,t / ) = G^(t-t',z) = B Y^ 



(-( t -t'^ie) 2 + z 2 ) 



- c e(t - 1') 



(2) 



+,- 
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Figure 1: The scalar field in the vicinity the segments labelled by =F contains divergence 
of ingoing/outgoing type. The segments of the trajectory are also labelled by re if they are 
connected to the boundary point t n . See the text for more details. 

where 

d+i i r(c)r(i) 

C = , -On = j — — . (3) 

2 ' 7T r(|)2 c v ; 

The presence of the mirror will scatter G R in a way consistent with its retarded nature. 
The net result we will find is that the divergences are propagated through the bouncing of 
G R between the mirror and the boundary. 

FigQ] illustrates the mechanism of the divergence matching method: A light ray 
starting at t' from the boundary follows a bouncing trajectory composed of segments of 
null geodesies. We label each segment by a number n and sign — /+. The number label n 
indicates that the segment is connected to the boundary point t n . Two segments joining at 
the same boundary point have the obvious interpretation as corresponding to the ingoing 
and outgoing waves, respectively, which we label by — and +. In the vicinity of each 
segment, G R is singular. In the following, we will simply refer to the divergences of G R 
as being of ingoing and outgoing type for the sign label — or +, respectively. The idea 
is to determine all the divergences along (n, +), (re, — ) through matching along the mirror 
trajectory and on the boundary. The matching procedure is explained in detail in the 
subsequent. 

The initial data for the divergence along segment (0, — ) is provided by G R , just as in 
[Hj. It is given by 

G R (t -t'^z) ^ ±B Q 

with c as in ([3]). Bearing in mind that the prescription t — > t — ie (t — > t + ie) corresponds 
to a positive (negative) frequency contribution, we may identify the upper/lower signs as 
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contributions from positive/negative frequencies. We may treat them separately in what 
follows. Using >/< for positive/negative frequencies, we write 

G °>> = {-t + t' + ie + z)^ G °>< = {-t + t'-ie + zY (5) 

The first step of matching is to be done in the vicinity of (to,zo = /(to))- Along segment 
(1, +), the outgoing type divergence of G R assumes the generic form 

G t> = fa-t-z + ie)*' G ^< = fa-t-l-iey (6) 

In the vicinity of (to;/(*o))> the divergence of G R is a superposition of Gq and G±. The 
Dirichlet boundary condition on z = f(t) gives 

Go,> + Gt,>=Q> C - < +G+ < =0. (7) 

Denoting t = to(l + x), then z = /(to) + f'(to)tox. Plugging J5J) and (JBJ) into (JTJ) and 
expanding to leading order in x, we obtain 

Bl =- B ji±mY- is) 



.W(fo). 

The next step is to do the matching at the boundary close to t\, where G R is a superposition 
of Gf and G7 . This will allow us to determine G7 . It is convenient to use the frequency 
representation. We consider the spatially integrated correlator which corresponds to the 
zero momentum mode. Therefore, for a single component of frequency we have, assuming 
ui > 0, 

GZ : zU-^H^iuz) = z¥,/I e --WKf) , w _> oo , (9) 

2 V VTW 



G- : zle-^H^i-cJz) = z^r A —e-^+^^-r) , w ->• oo, (10) 

2 V — VTW 



G+: zie~ iMt Hf{uz) = z^\ — e — , w oo , (11) 

2 V VTW 



G+ : z^e-^ff^C-wz) = A /-?- e -^-*(«*+¥) , w -> oo. (12) 



In order to make the boundary free of source away from t' , we need to require the super- 
position of Gi and G7 contains z? Jd(\uj\z) only, thus for each H^\cuz) in G>, we need 



2 



2 



H^\ujz) in G>. Similarly, for each Hy(-uz) in G <5 we need H^\—uiz) in G<. Using 

2 2 2 

the representation ([9]), we conclude that 

Bie~ ilTC _ -Bie inc 

G h>= fa-t + z + ie)*' Gl < <= {h-t + z-iey (13) 
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Repeating this process, we find the divergences propagated along the null geodesies are 
given by 



n '> {in-t-z + ief n '< (t n -t-z-iey 



G - > = 7= M? > G- n< = ~ Bn{e ] , , (14) 

n,> (t n -t- z + ie) c n,< (t„-t- z-ie) c ' 

where B n satisfies the recursion relation B n = —B n -\ ( i*/'^"""*) ) • Finally, we wish to 
extract the divergences of the unequal-time correlator again by taking advantage of the 
frequency representation, for which we also need the coefficients of the z d term in the 
expansion near z = 0, 



— jivd j 

> f V 7 isin^r(g + l) \2) y 1 

G< : zh-^H^i-uz) -> / f Z^V e --V, (16) 

< f v ; i Sin ^r(| + i) V 2 y v ; 

iird 



d 

i sm 



G+ : z^rtfH - J^^J— Q~ 2 e~^, (17) 



2 " V2 

ind , . d 



zsinf r(| + l) V 2 

Gn,>j Cn,o and G~< lead to four contributions to the divergence of G R (t — > t n ,t'). 

We calculate them one by one, 

B n e 



(t n — * + ^ + « e ) Jo V 

where g>{oo) is the weight for frequency w. Thanks to the ie prescription, the above can be 
identified as Laplace transform, which can be inverted to give 

pmjj r> —imrcAuitn, c— 1 

9 > M = — — m • (20) 

The contribution to the divergence of G R (t — > i n ,t') is then given by 

L du9>H Z^4v(Li) (I) 



d 

2 e -iw(t-ie) 



isinf r(f + l) 

f5 ne -¥ e -^ c i /i\i r(d+i) 



r(c) isinf^ r(f + l)V 2 / (-it n + it + e) 2c ' 



(21) 
(22) 



Similarly, 
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which implies 



_ ^ 7lu) R /3*717TC Jiwtn . ,C— 1 

- ^ no ■ < 24 > 



The contribution to the divergence of G R (t — > t n ,t') is then given by 



^ 9 <(^) . e "L./.., (^ ; e "" + "' 



'isin^r(f + 1) V2 
^e^e™ 710 1 r(d+l) 



r(c) isin^ r(^ + l) \2J (ii n -it + e 



\2c' 



(25) 



The calculations of the contributions from G^ > and exactly parallel those for 



> 



and G„ <• We obtain 



( vf^^e— i /iy r(d + i) 

r( c ) isinf r(| + i)V2y H*„ + it + e) 2c ' 1 ' 



iird 



lB ne - — e in ™ 1 fiy T{d+l 



n>< " T(c) isinf r(f + l)W (^-^ + e) 2c ' 1 ' 

Summing over the four contributions, we obtain a neat result for the most singular part of 
G R (t,t') as t -> t n : 



Gn{t ~* tn ' " 2 f r( c) r(| + 1) I H + «n + -) 2c " H + £ - -) 2c 1 ' ' 



Comparing (|28|) with the result obtained using the integral form of the divergence matching 
method in p3], we find perfect agreement. 



3 Divergence matching for a gravitational collapse model 
3.1 A gravitational collapsing shell 

In this section, we study a gravitational collapse model with the purpose of gaining further 
insight into thermalization of the dual gauge theory. The model has been described in 
detail in [11] in the quasi-static approximation. We recall the key ingredients here: The 
model contains a homogeneous shell collapsing under its own gravity. The shell separates 
the spacetime into the parts above and below. By above and below, we refer to the region 
between AdS boundary and the shell and the region between the shell and the AdS interior, 
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respectively. The corresponding metrics are given by AdSs-Schwarzschild and by pure AdS, 

„ -fdt 2 + dx 2 + dz 2 /f 
above : ds 2 = *- = — (29) 

, 2 -eft 2 + cK 2 + dz 2 

below : as = = , (oU) 

z z 

with / = 1 - Z - . (31) 

4 

The z/j is the position of horizon, which defines a temperature z^ = -iL. As we will see soon, 
in Schwarzschild coordinates, the horizon is always "protected" by the shell, but in Kruskal 
coordinates, the shell will be able to cross the horizon. Note that we have used tf for the 
time coordinate above the shell to distinguish it from the coordinate t below. However we 
choose the radial coordinate z to be continuous across the shell. We choose to parameterize 
the hypersurface £ traced out by the shell by r, x such that the induced metric on £ is 
given by 



, (It 1 + dx 2 
dSs = z(t) 2 • (32) 



The trajectory of the shell as measured by the coordinates above the shell is £/(t), z(t). 
Similarly below the shell it is described by t(r), z(t). The velocity of the shell is defined as 



U ~ dr ■ 



II 



(tf,0,z) above, 



u^ = (t,0,i) below. (33) 



Comparing (J29J) and (|32p . we obtain the following relations: 

fi f 2 - Z j = l,i 2 -z 2 = l. (34) 

The unit vector normal to the hypersurface £ satisfies n-u = 0, n 2 = 1. It can be obtained 
easily as 

= (—,6,ifzf) above, 
J J 

n" = (zz,6,iz) below. (35) 

The falling trajectory of the shell is determined by the Israel junction conditions 
|23j . We have used implicitly the continuity of the induced the metric in ()32j) . A further 
condition is given by 

[Kij - jijK] = K S i3 , {Kij}S ij = 0. (36) 
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We use Greek letters for spacetime coordinates and Latin letters for coordinates on S. 7^ 
is simply the induced metric and -fQ,- = n a I + r^Srsj ) is the extrinsic curvature. 



The square (curly) bracket denotes the difference (sum) of the quantities above and below 
the shell. Sij is the stress tensor of the shell. To proceed, we use the ideal fluid type stress 
tensor: = {e{z)+p{z))uiUj+p{z)^ij. Note that we allow the energy density and pressure 
to depend on the radial coordinate. 

It is straightforward to calculate the nonvanishing components of the extrinsic cur- 
vature above the shell as 

A T t — ; 9 , J^xlxl — J^x2x2 — J^x3x3 — o • \o<) 

zdz Z L Z 

The counterpart below the shell is readily obtained by taking / = 1. Applying (136p . we 

obtain the following constraint on the energy density and pressure, 

ik - ^ <*> 

We choose the equation of state to be e = ^, which gives e ~ 2; 3(o+i)_ inserting this into 
leads to the relations 



Vl+^ - vTTi^ ~ z 3 ^ , (39) 



Vl + i 2 + \[JT^ ~ z 1 - 30 . (40) 
To be specific, we fix o = |, which corresponds to a conformal equation of state, thus 



+ ¥ - yf]T^ = bz\ (41) 
with b a parameter. It is a short exercise to show that 




1 



bz 4 



if = = *L_. (42) 

With (|42|) . we can easily solve for the trajectory of the falling shell in the coordinate above 
the shell. A convenient way to parameterize the trajectory is to use the initial radial position 
z s , at which the shell starts falling with vanishing velocity, i.e. z = 0. This gives a relation 
among three parameters b, z s and Zh, 

bzj + -L = 2 . (43) 

This has the physical meaning that the temperature Zh is determined by the energy density 
b and the intrinsic scale z s . From now on, we will set z^ = 1, with the understanding that 
all other quantities are measured in units of the temperature. Since z s < % = 1, we may 
easily convince ourselves that | < b < 1. 
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Figure 2: The regions of parameter space marked by the number bouncings of the light ray 
on the shell. A smaller t' and z s always lead to more bouncings. We observe a good linear 
relation between — and t' along the boundaries of regions T n . 

Zs 

3.1.1 Bouncing light rays in the gravitational background 

It is useful to study the behavior of light rays in the gravitational collapse background 
before we compute the singular part of the correlator in the next section. The massive shell 
always falls with a speed less than that of light, thus naively the light ray will bounce forever 
between the shell and the boundary. This is however not true because of the warping, which 
freezes both the shell and the light at the horizon. For a given parameter z s and the initial 
time t' that the light leaves the boundary, there are only a finite number of bouncings 
before the light ray asymptotes to the shell trajectory. The problem of finding the bouncing 
trajectory of a light ray in the gravitational collapse background can be studied with simple 
numerics. Fig|2]shows a chart of the number of bouncings in the parameter space of intrinsic 
scale z s and the time parameter t' . We have not mapped out the regions with n > 2. The 
reason will be clear later. For a given z s , we may define T n (z s ) as the critical value of t' , 
beyond which an n-th bouncing does not occur. The T n are curves separating regions with 
different number of bouncings. In general, smaller t' and larger — lead to more bouncings. 

Zs 

Physically, a small t' means that the correlator probes an earlier stage of the thermalization 
process, while a larger — corresponds to a higher intrinsic scale the thermalization process 
starts with. As argued in [TT], the scale should be related to the saturation scale of the 
nucleus in the context of heavy ion collisions. 

One may wonder if the above picture is an artefact of the Schwarzschild coordinate. 
In fact, the picture is more transparent in Kruskal coordinates. Working with Kruskal 
coordinates defined by 





(44) 
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the AdSs Schwarzschild metric becomes 



ds 2 



-4r» 



Az 2 



-dt\ + dx\) + 



,2 ' 



(45) 



In terms of the Kruskal coordinates, the trajectory of the shell is determined from (|42p as 



t K 



x K 



bz I xk 



\-bz" 
b 




(46) 



In Kruskal coordinates, the warp factor cancels when the shell reaches the falling velocity 

dx K (| " bz 4 ) t K ~ J{\+bz±) 2 -AX K 



dt 



A 



bz 4 ) x K -J(l + bz*y-4t K 



(47) 



oo. The shell will cross 



It starts from initially at z = z s and asymptotes to 1 as z 
the horizon at finite Kruskal time. Naively, at z = 1, ^ — bz 4 = \J (\ + bz A ) 2 — 4, we may 
conclude that = — 1 according to (}4"7]h This is not true, however, because we also 
have tK = xk- As a result (|4T|) becomes undetermined, L'Hopital rule is needed in the 
evaluation. A more careful analysis shows a positive velocity when the shell crosses the 
horizon. Numerical integration of (j47j) also confirms that it is a monotonous increasing 
function of z. 

The light ray trajectory in Kruskal coordinates becomes trivial. It is simply given by 
^7 = ±1, therefore we expect an infinite number of bouncings between the shell and the 
boundary. For a given z s , the critical value T n {z s ) corresponds to the situation that the n-th 
bouncing on the shell occurs precisely when the shell crosses the horizon. For t' > T n (z s ), 
the n-th bouncing occurs below the horizon, thus the reflected will not be able to escape 
from the point of view of the Schwarzschild coordinates. 

3.1.2 Divergence matching method for a probe scalar 

We now turn to the spatially integrated retarded correlator 

G R (t,t') = J "d d ^xe(t -t')[O(t,x),O(t',0)] (48) 

in the gravitational collapse background. We choose to study the behavior of a probe scalar 
in this background and extract the correlator of the dual operator 0{x). We will first derive 
the matching condition of the scalar on the shell, and then apply the method of divergence 
matching to determine the most singular part of the correlator without solving the scalar 
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wave equation in the bulk. For the case of a scalar in the gravitational collapse background, 
the boundary condition on the hypersurface of the shell is more complicated, which requires 
a generalization of the divergence matching method formulated in the previous section. We 
will still work in general dimension d , which makes the structure more transparent. Finally 
we will set d = 4. 

The behavior of a scalar in the gravitational collapse background has been studied in 
[9] . The matching condition is simply the continuity of the scalar itself and its flux. In our 
particular example, it is given by 

f +' = * . (49) 

nf ■ V(j)f = n ■ V^> 

The first line of (|49f) implies uj ■ V0/ = u ■ V^>, which combined with (|49p gives the explicit 
relations 

ifd tf (/>f + zd z 4>f = id t 4> + zd z 4> 

fa f <l>f + iffd s <l>f = zdt<i> + id x <l> ■ (50) 

4>f = 4> 

Since we are interested in the retarded correlator, the scalar wave contains only the in- 
going component below the shell. The wave in frequency representation is z~z H A (ojz) or 



d 



(2) 

z~zH\ (—ujz) for positive and negative frequency modes, respectively. Working in the UV 



2 



limit as in |14| , we may use the asymptotic expansion of the Hankel function to show 

dt<P + = ^±<f>. (51) 
2z 

Plugging ()5ip into (|50p . we obtain a matching condition involving (f)f only, 

(*/ + jj d t f h + (* + ifWztf = + ( 52 ) 



where tf = — j — and t = yl + z 2 . All the quantities are to be evaluated on the shell. 
3.2 Method of divergence matching 

Now we derive a recursion relation for the most singular part of the correlator. We again 
refer to Fig(TJ As in the Dirichlet problem in Section. 2, close to each segment of the trajec- 
tory, the singular part of the scalar wave is either ingoing (for the minus sign) or outgoing 
(for the positive sign). Furthermore, working in the UV limit allows us to approximate the 
ingoing/outgoing waves by WKB solutions. From now on, we will suppress the subscript 
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/ when the region below the shell does not enter the discussion. We construct the WKB 
solutions 

fz dj 



which have the obvious property 



^(ujzVe^^TPJe-^, (53) 



(±d t + /9 z )0 T = ^0 T /. (54) 

We split the wave into ingoing and outgoing components (f)f = </>_ + Inserting this into 
([52]) and using ([54"|h we obtain 



/ (VTT¥-^fT¥ ) d _ 1 

2d t J> + = — ^ : ^ — (0_ + 0+) . (55) 

V/ + z^ + z ^ 

This is the key equation for divergence matching. 

We begin with the matching at (to,zo). As we argued before, close to the segment 

0, — , the singular part of the wave is ingoing </>o,-> and close to 1,+, the singular part 

of the wave is outgoing 4>i t +. There is an overlap region between the two, which is close 

to (to,zo) on the shell. Since 4>o- is simply the WKB limit of thermal bulk-to-boundary 

correlator with ingoing component only, a WKB computation or simple analogue with zero 

temperature bulk-to-boundary correlator gives 



d-l 

z 2 B 



1 1 



-t + t' + ie + f^Y (-t + t'-ie + j, 



z dz' 



o 7RT 



%,- + <I>q,-, (56) 



with 



d-l d-l 

z 2 Bo < z 2 Bq 



-t + t' + ie + ' " f-t + t'-ie + J* 



(57) 



10 f(z')J V J ° /(*'), 

The contributions with a > and < sign have a clear interpretation according to the ie 
prescription: They correspond to the contributions from positive and negative frequency, 
respectively. Since we are interested in the most singular part of we may drop the 

term proportional to 0x,+ on the right hand side of ([55]) . It is easy to guess the general 
form of the singular part to be 



d-l 

z~Bi 



-t + t' + ie + f Q -fa 



(-t + t'-iz + fo W) 
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Inserting the above into ([53]) , and substituting t = to(l + x), we obtain for x — > that 



/(yr+i*-v7 + i 2 ) d _i Bo 2Bi(c-l) 



77+i2 + i 2z x x ' ' 



-x H , 2 x db e —x , z x ± e 



which gives 



Bl - B »4^ri) — iTFTI^lJ ' (60) 

The second step of matching is to be done near the boundary point t = i%, which is 
the overlap region of the segments (1, +) and (1, — ). The form of the most singular part of 
_ is determined such that <p = <f>i + + <fii _ is free of source, i.e. (f) contains only the vev 
term z d there. It is useful to separate the positive and negative contributions and look at 
them in frequency representation. The near-boundary solution of <j>± are given by 



For uj > : 



z$hP(uz), <t> + ~z*H®( U z), (61) 



For uj < : 



~ z*Ht\-uz) , </>+ ~ ^^(-u,*). (62) 



2 



In order to cancel the z° term near the boundary for the positive frequency contribution, 
we need to have a H^\ujz) in (j>\ — for each H^\ujz) in </>i,+- On the other hand, the 

2 2 

coordinate representation 4>i + and take the form 

>f + = ^ T , = ^ r> (63) 

' -t + t' + ie-S^) C (-t + t> + ie + $* ^y 

with # to be determined by the matching. Note that Jq j^jj z as z -> 0. Using the 
asymptotic expansion of the Hankel function, we conclude that 



("* + t' + ie + f* 



7RI 



(64) 



A similar analysis of the negative frequency contribution leads to 



d-l 



(65) 
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Repeating this process, we obtain for the most singular part of the correlator 



-t + t n + ie - J Q joy ] ( — t + t n — ie — J jAv 



f(z')J V " J0 f( z ') 

z^e- i7TCn B n t< z^e i7rcn B r , 



(66) 



"■ — / _ z d ' \ c ~"' ^ n > — / _ z d ' N C_ n ' 

i + i n + ?e + J j^tj J pi + in - «e + J Q y^y 

£> n is to be determined from the recursion relation 

where all the quantities are evaluated at t = t n -i, z = z n -\. 

Now that we have obtained the most singular part of <p n ,± close to the segments 
(n, ±), the next step to reconstruct the most singular part of the boundary retarded two- 
point correlator in the vicinity of t n . This is done as follows: We first write (p^± i n frequency 
representation, 

u > : 

f'°° d 

<f>> + = / dug> + (u;)z2H { ?\oJz)e- 1 " 1 , 

JO 2 

</>>_ = / do;g> t _(cj)z^H^\ U z)e-^ \ 

JO 2 

uj < : 
V+ ~ 



OO 2 



■ - f du J g<(u)ziH^(-uz)e-^ t . (69) 



OO 



In the UV limit, we can use asymptotic expansion of the Hankel functions on the right 
hand sides of (|68p . which become Laplace transforms (for negative frequency contribution, 
a change of variable is needed). These are easily inverted to give 



Sj+H = s£-M = -W-^e^e'^e-fe J " r ""' , (70) 

V 2 r(c — n) 

with c given by ([3]). Finally, the most singular part of the two-point correlator is obtained 
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by taking the coefficients of the z d term and convoluting with the weight in (|70p , 



G> + = [ du - ( 1 e— £ + ( w ) 

• /\,. fHP 1 /-V. 



Gl-= I du— ___(-) 5 >_H, (71) 



o 



o 



iivd 

e 2 


1 


isin ^ 


r(f + i) 




i 


—e 2 


• „• 7rd 

z sm — 






i 


— e 2 


isin ^ 


r(f + i) 






e 2 


i 


isin ^ 


r(| + i) 



-oo 



= / ^ ^7-TT77^ „ ( ) e-^^< ( W ) . (72) 



It is easy to perform the integral and we combine the results from positive and negative 
frequency contribution as follows, 

r> = r> +r> r(2 C -n) B n e~™ ^ 

r< = r< + r< ^ r ( 2c ~^) V^ 11 m s 

" n ' + 2-ir(f + 1) r(c-n) (-t + t„ + ie) 2c -"' 1 ; 

Note that the singularities of the correlator become milder as n increases. In other words, 
every bouncing of the light ray on the shell lower the power of the singularities by one. 
This is a consequence of the same effect on the bulk scalar close to the segments depicted in 
FigJTJ Furthermore, there is also a suppression on the numerical factor B n by each bouncing 
as dictated by §6j 



3.3 An explicit example: The quasi-static state 

In this section, we will test our algorithm with an example where explicit analytic compu- 
tation is possible. Since the procedure outlined in the previous section is valid for all shell 
trajectory, we choose to work in the case of static shell: z = 0. This can be viewed as a 
quasi-static state studied in [TT] . 

Since the shell is static, modes with different frequency modes decouple, which allows 
for a simple treatment in frequency representation. We start with the matching condition 
(f49l) . Setting z = 0, we have 

H = <t> | h = yjt (74) 

yffd z 4> f = d z (f> I Vfdz(f>f = ^j d z4> 

The right hand side is the Fourier transform of the left hand side. We have used a tilde 
to indicate quantities in the frequency space. When i = 0, the mismatch between the 
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time coordinates above and below the shell is simply: dtj = -yrdt. This also leads to a 
corresponding mismatch in frequencies and the additional factor -^j on the right hand side 
of (jTSj) [11] . We will take the frequency above the shell to be u, thus the corresponding 
frequency below the shell is -^j. 

In the UV limit, <fi satisfies d z <p — ^jcj) = ^-4>- Plugging this into ([73]) . we obtain 

fdz4>f ~ iuh = ( 75 ) 
We again split the wave above the shell into ingoing and outgoing components, 

0/ = + B4> + , 

where (j) T = z^e^ ^ H^ 7 ) . (76) 

The ratio is readily determined from (|75p . To further simplify the calculation, we keep 
only terms to the leading order in / (or 1 — z s , where z s is the radial position of the shell) . 
We obtain 

A -Aiuj </>+ -Moo -2iu>$*J^ 



b v7(rf-i)0- v7(rf-i; 



e — Jo 7(F) _ ( 77 ) 



To calculate the two-point function, we need an expansion near z = 0. Note that (f>± defined 
in (|76p is not valid as z — > 0. However, they may be matched to the near-boundary solution 
which is given by Hankel functions, 



For uj > : 

: -r(>) • • . -ln( 2 )< 

d 

2 2 

For oj < : 



Z 2H); ) (ujz)e—, (f> + ^ Z2H { d 2) (uz)e— , (78) 



-iwc d _ j^^J s iirc 



z2H y d >(-ujz)e—, (f) + ^ z2H\ l {-ujz)e— . (79) 



2 



With (|77p and the mapping above, we may write down the retarded correlator for separate 
contributions from positive and negative frequencies as follows, 



■ 71 Li , '171 (J ■ 71 Li 171 Li. _, / ^ ,J \ 

> _ -L f^-^ -^-^^^r f „ r _ (go) 



; TTtZ J ITTC „■ 7Tt2 ^ — ItTC / „ ^ 

2^ 7o Ae^ - Se'T T(l +|) V 2 

A direct integration of ([80]) is not possible. Note that 4 is suppressed by -j- in the UV limit. 
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We may perform a series expansion in of the integrand using 



2tt 



G> = J- / due-^'U-e-^ 



_i7r<2 



2vr J_oo 



n>l 



n>l 



5 

I 6 " 



r(l-|) /a; 

r(i + g)b; ' 



(82) 



-id 



r(i-|) 
r(i + |) V 2 



(83) 



The lowest order terms, i.e. the 1 in the square brackets, simply give the vacuum retarded 
correlator. The n-th order terms precisely give us the most singular part of the retarded 
correlator, 

T(d - n + i) e -iTc(n-i) 



isinf T(l-f) fVf(d-l) 

^nfi — : ~ 1 " 



2 d vr r(l + f) V 4 
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d-n+1 ' 



zsinf r(i-f) / V7(rf-i) ' 



T(d - n + i) e ^(^-i) 

d-n+l 



if 



(84) 



'o 7RT 

We note that the dependence on the temperature is not visible for the WKB approximation 
we are working in. 

Let us calculate the same quantities by applying our algorithm outlined in the previ- 
ous section. The only input required for (jT3|) is the B n . Working to the lowest order in /, 
we obtain from (|67|) 

d-1 



B n 



-B n -iy/f . 



4(c - n) 

The initial condition for the recursion equation is Bq 

dz' 



(85) 



- , d \ 2 . Noting that t T 



2 c 7r r(f) 



t' + 2n Jq jjjpj in case of a static shell, we can show that (I84p and (I73p are in perfect 
agreement! 



3.4 Thermalizing state 

Now that our algorithm has survived the non-trivial test, we can apply it to the more 
complicated case of gravitational collapse, with trajectory of the shell satisfying 



dz 
dF f 



1 



i>~ i. 



(86) 



We may deduce the properties of the most singular part of the correlator already from (|73|) . 
We will specialize to d = 4 from now on. For a given time t' , the spatial integrated retarded 
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correlator / d 3 x 9(t - t')[0(t'), 0(t)] contains singularities as t — )■ t n , with t n being the n-th 
boundary point on the trajectory of the bouncing light ray. The most singular part of the 
correlator takes the generic form 

A — I7rc(n— 1) a iirc(n-l) 

W (-t + in + ief-^ n ~ (-t + in + ief-n { ' 

The numerical factor A n has to be determined for a given trajectory. The condition of the 
divergence of (j66|) implies n < c = |, which is the reason we have mapped out the number 
of bouncings only upto n = 2. The n = case gives simply the lightcone singularity. The 
nontrivial singularities appear for n > 1. We have shown in Section 3.1 that i n — > +oo as 
t' — > T n (z s ) from below, which means the retarded correlator will be free of singularities 
when t' > T\{z s ). Therefore we propose to use T\{z s ) as our definition of the thermalization 
time ith- Let us discuss the connection to heavy ion physics. Restoring units, we have 

tth = X \ T S \ (88) 

where T corresponds to the initial temperature of the quark-gluon plasma (QGP) formed 
in heavy ion collisions and z s is related to the saturation scale of the nucleus. We see 
that the thermalization time is inversely proportional to the initial temperature and the 
numerical factor T\ involves an interplay between the temperature and saturation scale, 
which is ultimately determined by the collision energy and size of the nucleus. 

In Kruskal coordinates, the light ray will always have at least one bouncing. But 
for large t', the bouncing will occur only after the shell crosses the horizon, thus the light 
ray will not be able to return to the boundary. From the point of view of Schwarzschild 
coordinates, this corresponds to the case where the light ray never reaches the falling shell. 
Therefore for sufficiently large t' , we expect the singularities of the correlator to be identical 
to the singularities of thermal correlator evaluated in the background of an eternal black 
hole. In the next section, we will find the singularities in thermal correlator actually appears 
in the complex t plane. 



4 Singularities in thermal correlator 

We now move to the second part of this paper where we establish a relation between the 
singularities of the retarded correlator in coordinate space and the quasi-normal models 
(QNM). We first do this for a field theory in thermal equilibrium, which is dual to the 
AdS-Schwarzschild background. We obtain explicit expressions for thermal field correlator 
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Figure 3: Bouncing of a null geodesic in the AdS black hole in Kruskal coordinates. The 
geodesies starts on the AdS boundary on the right and bounces off the future singularity, 
the second AdS boundary and the past singularity back to the boundary. A similar set-up 
has first been discussed in |22j . 

in coordinate space. These expressions allow us to locate singularities in the complex 
time plane, which turns out to be closely related to the QNM. The first indication of this 
connection was discussed in [22], where a geometric optics picture has been used to derive 
the asymptotic QNM. It was assumed that the coordinate space correlator has periodic 
singularities in complex time, with a complex period determined by the time a light ray 
needs to make a loop by bouncing in the full Penrose diagram, as shown in figure [3l 

Applying this picture here, we may draw the following conclusion from our study 
of the gravitational collapse model: As we increase the initial time t', the unequal-time 
correlator probes different stages of the thermalization process dual to the gravitational 
collapse model. The resulting singularities of the correlator show a shift from the real time 
to complex time. We will solidify this picture by explicit evaluations of the correlators in 
the BTZ black hole and AdSs-Schwarzschild black hole. We will see that singularities are 
indeed present in both cases. 
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4.1 Unequal time correlator from BTZ black hole 

We are interested in the singularities of correlator in coordinate space. The black hole 
background preserves time translational symmetry, therefore the unequal time correlator 
only depends on the difference of the insertion times of the operator, G R (t, t') = G R (t — t'). 
Below we will set t' = for convenience. We start with the correlator in momentum 
representation, 

G r (uj) = J dte iult G R {t). (89) 

We will obtain G r {oj) by solving massless scalar wave equation in the BTZ background. In 
the reduced unit 2ixT = 1, the BTZ black hole takes the form 

z 2 

where <j) is periodic eft ~ (j) + 2tt. The solution of the massless scalar satisfying the ingoing 
boundary condition at the horizon and approaching unity at the BTZ boundary is expressed 
in terms of a hyper geometric function, 

$( Z ) = -1 2V H_ z 2\ 2 z 2p I l l ;l-iuj;l- Z 2 . (91) 

v ; T(l - iuj) v ; V 2 2 / v 7 

Note that we consider a solution without ^-dependence, which corresponds to the spatially 
integrated correlator where an integration over 4> is implicit. ([9T]) leads to the retarded 
correlator 

G R (u) = iuj - - y ). (92) 

All the poles of Gr(uj) lie in the lower half plane at w = — i(2Z + 2), which means 
the correlator is nonvanishing only for Ret > 0. In the latter case, we close the contour 
clockwise and pick up the poles in the lower half plane. It is easy to sum over the residues. 
We obtain 

G R (t) = -2 Y e-^\2l + 2f = (93) 
fo sinh 1 

Prom the picture of a light ray bouncing in a confining box of the Penrose diagram 

as shown in figure [3l we expect singularities to have a period of At = itr, which is the time 

a light ray needs to finish a loop. Explicit expression for the correlator (|93p indeed shows 

singularities at t = mn for n £ Z. Strictly speaking, we need to give t an infinitesimal 

real part to justify Ret > 0. This will not be necessary for the case of AdSs-Schwarzschild 

black hole below. This difference arising from the global structure of spacetime between 

BTZ space and higher-dimensional AdS-Schwarzschild space has been discussed at length 

in [25]. 
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4.2 Unequal time correlator from AdSs Schwarzschild black hole 

Now let us look at the correlator in the AdSs Schwarzschild background. The unequal-time 
correlator is again given by 

G R {t) = -!- [doj e" 1 ^ G r (uj). (94) 



2vr „ 

The derivation of the relevant part of G r {oj) can be found in the appendix. Again it is 
nonvanishing for Ret > 0. We close the contour clockwise and pick up poles in the lower 
half plane. The poles from the QNM is expected to give singularity of G R (t) in the complex 
t plane. Since the explicit expression for G R (u) is not available, we will try to evaluate (|94p 
approximately. This amounts to summing over the contributions from poles with large \co\. 
The analytic expressions of asymptotic QNM for AdSs Schwarzschild have been obtained 
in [261 [28]: 

nir + 6 nil + 6 , . 

u n = , co n = — : , (95) 

Xq Xq 



where xq = xq = i and 6 = ^ + ^^PH To proceed further, we need to know the 
residues at the QNM, all of which are simple poles. This can be achieved by generalizing 
the work [26J. We collect necessary steps in the appendix. We obtain for the residues of 

mK = =! + ») = , res(k ,„ = =!L + » ) = -i* (96) 

Xq OZXq Xq OZXq 

Combined with the fact that G r {uj 3> T) ~ a; 4 lncj, we make an educated guess for G r (oj), 
G>) = 4^(^) + ^=)V (97) 

61 \ 7T 7T J 

This approximate expression is also consistent with the symmetry Gr{<jj) = G* R (—u*), 
which is present for \u\ ^> 1, when we may ignore the dependence on the offset 8. Now we 
are in a position to evaluate (194p with the approximate expression (|97p . 



The contribution to G R (t) from the QNM at calculated to be 

i sr^ i, , i a 7T i _i(™+B)t / n7T _|_ 0\ 4 n 
32^ n x 32 V / xq 

= L e ~% (ly^e-n^t) A^-e^Oi \ (98) 

32 V^o / k 32 {it^2m x XQ%f' v ; 

where $ is the Lerch transcendental function and rri\ £ Z. The symbol — means that 
the most singular parts of the two sides are equal. The contribution from the QNM at 



#1 Our definition of xq differs from that of [26] by complex conjugation 



23 



to G R (t) can be easily obtained by the substitution xq — > xq in ([98]) . with the 



singular part 




n 



Xq 32 



(it + 2m2Xoi) 



1 



5 - 



(99) 



where m2 £ In the region Ret > 0, we find singularities at t = —2miXQ and at 
t = —2m,2Xo for m\ < 0, m,2 > 0. The singularities indeed have periods dictated by the 
dual periods of the QNM ^- and J^-. This completes our discussion of the singularities 
in thermal field correlator and confirms the existence of singularities in the complex time 
plane. 

5 Singularities in unequal-time correlator from evolution of 



The explicit calculation in the previous section has taught us an important lesson: The 
singularities of the correlator in the complex time plane are closely related to the structure 
of the quasi-normal modes (QNM). In this section, we aim at improving our understanding 
of the singularities of the unequal time correlator by studying the evolution of QNM in 
the gravitational collapse process. This is only possible for the quasi-static state, which is 
defined by a sequence of states corresponding to the shell held at different positions above 
the horizon. As the shell is lowered to the horizon, we expect to observe features of the 
thermal field. While the sequence of states does not correspond to physical state undergoing 
thermalization, it does help us to gain more insight to the thermalization process. We 
specialize to d = 2, corresponding to the formation of BTZ black hole, which we have a 
better analytic control on. We study the evolution of the singularities in the unequal time 
correlator from the dual evolution of the QNM. We reproduce the singularities of (|84p in real 
t from the normal modes, which are a subset of the QNM. We also explore the possibility 
of having singularities in complex t from the other QNM. 

For the quasi-static state, we can look at the retarded correlator G q R in momentum 
space, with the superscript q indicating it corresponds to a quasi-static state. We will study 
its QNM in the whole complex plane. We need to work out the matching condition on the 
shell in momentum space, which is valid for arbitrary complex uj. The starting point is 



QNM 



again the RHS of (JHD, 




(100) 
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with to and -Hj being the frequency above and below the shell respectively. For complex u, 
the scalar below the shell is given by 

= ^i(e-f-^=z). (101) 
We may check that for oj > 0, this reduces to the real expression zH-jp (^jz). Conversely, 

(2) 

for uj < we find zH\ (—-^=z). There is a branch cut on the negative imaginary axis, 
which is crucial in reproducing the retarded correlator in coordinate space. 

The scalar above the shell is a linear combination of ingoing and outgoing waves 
(pf = Acj)_ + B(j} + . When shell is sufficient close to the horizon, ^>=p is approximated by a 
series solution near the horizon. We choose 

4> T = (l-z)^e ±iujlJ ^ +0{l-z). (102) 

The z-independent factor is introduced such that (ft T reduces to (|76p as z — > 1. Now inserting 
this into (jlOOp and performing a series expansion in / or effectively in 1 — z, we find that 
the corrections to series solution cj)^ are subleading and an asymptotic expansion of (jlOip 
gives 

i 

l(jJ 1 I !W 1 / ITT \ O , 3 



4 



J B = eV7/- 1 -^ i ^ . 103) 

We note that although the branch cut along the negative imaginary axis exists in both A 
and B, it is absent in their ratio 

The ratio agrees with what we obtained for real oj in (1770 in the limit / — > 0. A more 
general expression can be obtained for an arbitrary choice of the ingoing/outgoing wave 
solutions in d dimensions, 

b = Hd-DVTi. 

A 4w + ' 

to be evaluated at a position of the shell sufficient close to the horizon. W only need (|104p 
for our purpose, from which we may derive the corresponding retarded correlator for a 
quasi-static state, 
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where the superscripts and d denote the source and vev of the corresponding quantities. By 
construction, we have (fi'L = Gr{oj)4P_ and <p^_ = Ga{<^)<P% up to contact terms. Therefore, 
we obtain 

G R {«) - • (107) 

Now we use exact ingoing/outgoing wave solutions in BTZ black hole background, 

h = (1 - ^ **F (l T f , 1 T y ! 1 T to; 1 - **) . (108) 

The normalizations in (|108|) are chosen to match the near-horizon behavior in (I102p . It is 
not difficult to obtain from (|108p that 

~ r (i — iuj) ~ r (i + iw) 



r(i-T) 2 ' + "r(i + f) 2 ' 



(109) 

= iw - w 2 V \1 ~ y J , G A (u) = -iu - ^ f 1 + y 1 , (110) 

where we have suppressed the overall numerical factor in Gr and G^- Inserting (|109|) and 
(fTMj) into (HUZD, we obtain 

G * {uj) = 2 r ( i-^) r(i+L) ■ ( 1U ) 

r(i-^) 2 ^ 4lw r(i+^7 

Now we are in the position to investigate the structure of QNM corresponding to the 
retarded correlator (|llip . First we note that possible QNM from Gr{uj) and G a (uj) at 
uj = 2i{n + 1) and oj = —2i(n + 1) do not arise since the prefactors r } 1 ~^}-2 and r ( 1+ta; ) 



r(i--f) r(i+f 
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evaluate to zero, giving rise to a finite result. Therefore, we conclude that the only QNM 
arise when the the denominator vanishes. It is useful to write (jllip as 



* _ r(i-f ) r(i+M 

r(i +¥ ) 2 r(i-M _ 
r(i-i^) 2 r(i+j W ) J 4luJ 



Let us consider asymptotic QNM with \ui\ S> 1. The limit lim a ,_ i>00 — -, — ^-4 depends 

11 r(i-2f ) r(i+iw) 

on the argument of w. The role of the dependence of the correlator on the argument of ui 
has been emphasized in [27]. We first look at the limit along the real axis, 



r(i + f) 2 r(i-i^) 



r(i-ff r(i + 



IUJ 



i4 luJ oj — > +oo 

(113) 

-i4 4tJ w -> -oo 
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Figure 4: The left plot shows 



r(i+f ) 2 r(i-^) 



r(i-2 



r(i+iw) 



versus and the right plot shows / luJ jJ- 



versus iu. 



Using ()113p . we find the asymptotic roots u = a + ib with \u\ 3> 1 for fixed /, 

(2n - 1)tt 



ln| 

hT|' ~ln~T~ 



hrf 
ln{ 



(114) 



with n large positive integers. We find in both cases - — > as n — )■ +oo for fixed /, thus 
the roots are asymptotically real, which justifies the use of (|113p . We identify the first 
and second lines of (|114p as positive and negative normal modes. These normal modes 
are equidistant with a spacing Au = ~~j~ri which gives rise to the period At = — In j in 
the singularities of the unequal time correlator in coordinate space. Note that the period 
is robust in that it is only governed by the phase factor ^\ . On the other hand, the 
period according to (|84"j) is given by 2 J* jfg^ — m (l — -z) + In 2 — ^ — In ^ as z — > 1. We 
see that the two approaches give consistent results. 



Next we look for QNM along the imaginary axis. For purely imaginary u, the de- 
nominator is guaranteed to be real. It is convenient to study its behavior by plotting 



r(i+^)-T(i-Mj) 
r(i-^) 2 r(i+iu/) 

and too = —2n and zeros at iui 
has the simple asymptotics 



and / as functions of iu. The former has infinities at iu = 2n — 1 

-(2n — 1) and iu = 2n with n positive integers. The latter 



1 Aiu 



iu 
+oo iu 



-oo 
-oo 



(115) 



Fig. U] shows them in separate plots. It is easy to see that the QNM are given by their 
crossing points at approximately iu = —{2n — 1) and iu = 2n — 1. The appearance of 
QNM in the upper half plane might seem worrisome at the first sight. In fact, this is 
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just a reflection of the fact that our quasi-static state is unstable against evolution toward 
equilibrium: From the bulk point of view, we need an external force to keep the shell 
levitating above the horizon. Note that for the thermal field, the QNM in the retarded and 
advanced correlators are at iu = 2n and iu = —2n. Our QNM of quasi-static state agree 
with the thermal QNM in the period. However, they differ in the absolute value by an 
offset. Note that we are studying the QNM in the retarded correlator for the quasi-static 
state. Nevertheless, we have found the QNM corresponding to the retarded correlator at 
approximately iu = — (2n — 1) and iui = 2n — 1, which are democratically distributed along 
the positive and negative imaginary axis, i.e. there is no sign of the retarded nature of the 
correlator. This is only visible when we study the residues at the QNM. Since all the QNM 
are first order poles, the corresponding residues are obtained by taking the derivative of the 
denominator with respect to u. The residues of G q R {u) are given by 



i res 



r(i-f ) v(i+iu) 



d(iui) 



r(i+f ) r(i-fcQ _ iu] 
r(i-^) 2 r(i+i^) J Aiu) 



r^jj (Gr(v) - GUM) 
'r(i+f) 2 r(i-M f _ iujVT 



(116) 



We have inserted an i to the residue for convenience. The residues are found by first 
obtaining the QNM by solving 

r(i + f) 2 r(i-M _ _v7 



r(i-f) 2 r(i + 



IU 



and inserting them into (|116|) . Table [T] shows the residues for several QNM, which encode 
the retarded nature of the correlator: The residue at iu ~ — (2n — 1) is always smaller than 
that at iu ~ (2n — 1) and their ratio decreases further as / — > 0. 

With the above analysis, we have shown the existence of QNM along the imaginary 
axis at symmetric locations, but with asymmetric residues. Moreover, we expect by in- 
specting the properties of QNM along the real axis that there are also singularities in the 
complex t plane with period At = in. However, it appears that in the WKB approach used 
here, the possible singularities in the complex time plane are missed out. Note that the 
WKB approximation in Sectionj3]3] is valid for real oj only. 

To see if singularities in the complex t plane indeed exist, we need to work a little 
harder to obtain the residues for the QNM. Analytic calculation is possible asymptotically 
for large QNM with n>l, due to the fact that the QNM asymptote to ioj = ±(2n — 1). 
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£=0.9 




n=2 




n=3 


n=4 


% ies(iu rj — (2n - 
i res(zw « 2n — 


-1)) 
1) 


-2.0 x 10" 
0.72 


-5 


-4.6 x 10~ 8 
9.3 x 10~ 3 


-1.1 x 10~ 10 
8.4 x 10~ 5 


ratio 




-2.8 x 10" 


-5 


-4.9 x 10~ 6 


-1.3 x 10~ 6 


f— n 1 

1 U. 1 




n — 9 






n — A 


i res(iu; ~ — (2n - 
? res(iw ~ 2n — 


-1)) 
1) 


-4.2 x 10" 
8.8 x 10- 


12 
6 


-1.5 x 10~ 18 
2.3 x lO" 11 


-5.4 x 10~ 25 
3.3 x 10~ 17 


ratio 




-4.8 x 10" 


-7 


-6.2 x 10~ 8 


-1.6 x 10~ 8 



Table 1: The residues of G q R {ui) at iui « ±(2n — 1) for several n at / = 0.9 and / = 0.1. 



At large n, the locations of the QNM are 

c2n- 



iu: = — (2n — 1) 



r 



4(2n - 1) r(| - n) 2 T(2n)(2n - 2)! ' 



ioj = (2n- 1) -4(2n-l)/ 



In- 



T(I + n) 2 



r(| - n) 2 T{2n){2n - 2)! 



(118) 



Note the last term in both expressions in (|118p are corrections, which tend to zero rapidly 
as n grows. Working to the lowest order in the corrections, we obtain the following results 
from (fTT6j) . 



For ioj ~ — (2n — 1) : 
f—iui v7 



r(i+i f ) r(i-fc,) _ ^ 
r(i-i^) 2 r(i+j W ) J 4luJ 



For iw ~ 2n — 1 : 

f—iu> v7 



a!(iu;) 



r(i+^) r(i-w) 
r(i-^) 2 r(i+^) 



r(i + n) 2 



4/i" 2n (2n - l)r(| - n) 2 r(2n)(2n - 2)! 



(119) 



4(2n - 1)- 



T(I+n) 2 



/2- 2ra r(| - n) 2 r(2n)(2n - 2)! 



(120) 



The remaining factor Gr(uj) — Ga(oS) in (|116p is also linear in the corrections, with the 
following expressions: 

^2- 



G R {u) - G A (oo) = (iu> ± (2n - 1)) ( -y ) (2n - l) 2 



7T 



(121) 
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Inserting (|119p . (|12ip and f|X18j) into (|116p . we obtain the residues 
For « — (2n — 1) : 

ires = 5 — ^ i , 122 

32 \ v r(|-n) 2 r(2n)(2n-2)!/ 

For iu) 2n — 1 : 

ire s = 8(2„-W ^!Z£li±^ V (123) 
We may approximate the factor involving gamma functions as follows, 

Km r ^ + "' 2 ~» (124) 

n ^+°°r(| -n) 2 T(2n)(2n-2)! ^ • 2 4 ™ 

For Ret < 0, we obtain the contribution from the QNM on the positive imaginary axis, 

G%(t) » V e-^iresfiw » -(2n - 1)) = — 7 — r-. (125) 

* V 2 7fe 2i+41n ^-l X 



For Ret > 0, the contribution from the QNM on the negative imaginary axis has an 
additional enhancement factor (2n — l) 4 , and we end up with 

G q R it) « e-^iresCiw « (2n - 1)) = ^^(e^ ({) ' , -4, i). (126) 

With our final expressions (|125p and (|126p . we may answer the question about the singular- 
ities at complex t: Possible singularities in ()125p and (|126p are at 2t + 41n ^ + 2nim = and 
at -2t + 41n{ + 2irim = 0, respectively. However, since In { < 0, both cases are excluded 
by the corresponding constraints Ret < and Ret > 0. We therefore conclude our search 
for the singularities with the conclusion that for the quasi-static state, the singularities in 
the complex t plane for are absent. 

We leave an analysis beyond the quasi-static approximation for further work. Our 
result for the quasi-static case is very different from a physical state undergoing thermal- 
ization, as described by the moving shell. The shell shows distinct trajectories in Kruskal 
coordinate for the quasi-static and moving shell cases: In the quasi-static state, the shell 
"moves" along x\ — t 2 K = constant, while the moving shell takes a trajectory ■^ s - = — 1. 
The difference is analogous to that of a free-falling object and a standing, but Unruh ac- 
celerating object [29J. It remains to be seen whether singularities emerge in the complex t 
plane for a physical thermalizing state. 
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6 Discussion 



We have generalized the divergence matching method for a gravitational collapse model. 
Application of this method leads to predictions on the evolution of the singularities in 
unequal time retarded correlator G R (t, t') = 6(t-t')[0(t),0(t')} for a state in thermalization 
process. The singularities of t are dictated by the geodesies of a bouncing light ray initially 
starting at t' . The fact that we found singularities at real t can be traced to the nature of 
the divergence matching method, which is essentially a WKB approximation for real ui. We 
have shown with an example of quasi-static state that possible singularities may exist also 
in the complex t plane. It will be very interesting to confirm this, possibly by extending 
the current WKB approximation to complex u). The singularities at real t coming from 
the normal modes have clear separation between contributions from positive and negative 
frequencies as shown in (|84p . This will not be the case for singularities at complex t coming 
from the QNM, in which no ie prescription is needed. This is related to the fact that in the 
positive and negative frequency in the black hole background is defined with respect to the 
Kruskal time. It will be interesting to see if the evolution of singularities at complex t can 
shed some light on this. 

Another interesting point to note is that the formulation of Sectionj5]does not require 
|w| S> 1. Actually it is sufficient to have ^ 3> 1, which is always realized when the shell is 
sufficient close to the horizon. However there is one exception: the gapless hydrodynamic 
mode. For the case of the dilaton, the QNM does not contain a gapless mode. However a 
gapless hydrodynamic mode does become relevant for example for a probe gauge or gravi- 
ton field, due to the presence of conserved charges |30| . It will be interesting to study the 
evolution of the hydrodynamic mode. We hope to return to these issues in the future. 

Note Added: When the paper is near completion, we received [31], which has partial overlap 
with Section 5 of our paper. Qualitative agreement has been found between the QNM that 
are asymptotically real (normal modes in our paper and blue points in their plot). The 
counterpart of their red point seems to be absent in our paper. The particular QNM may 
be related to the hydrodynamic mode unique in their model. 
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A Residue at the QNM 

We closely follow section 3.3 of [26] in the calculation of the residue. Readers are encouraged 
to refer to [26] for more details. The tortoise coordinate is defined through x = /J j^j^ 
with an appropriate choice of branch cut. Near the boundary r = oo, the solution to the 
wave equation is given by 



&(x) ~C + v / 27ru;(2 - xq)Jj_oo (oo(x - x )) + C-W2mJ(x - xp) J-j^ (u(x - x )) 

2 2 

~ (c+e ll3 + + C_e i/3 ") e *"G*-*o) + (c +e -^+ + CLe^M e -W*-*o) ^ (127) 
where f3± = \ (1 ± joo). Near the singularity r — 0, the solution takes the following form 

ttujxJ 3_{ojx) + B- V 2ttujxJ_ 1 (ujx) 

2 2 

~ (B + e~ ia + + B_e~ ia -) e iuJX + (B + e ia + + B_e ia ~) e~ iuJX , (128) 

where a± = ?(l±j). (|127p and (|128p can be matched along the same Stokes line. To 
impose the ingoing boundary condition at the horizon, which lies on another Stokes line, 
we need to rotate (|128p to the same Stokes line as the horizon. After the rotation, (1128P 
becomes 

~ (B +e - ia + + B_e~ ia -) e iuJX + (B +e - 3ia + + S_e" 3ia -) e" iaJx . (129) 
The ingoing boundary condition gives 

B +e - 3ia + + 5_ e - 3ia - = 0. (130) 
Matching (fTZTj) and (fT28|) and using (fT30|) . we obtain 

£f ^g2ia+ g2io_ ^ giajx— i/3-f _|_ ^g4ia+ g4ia— \ g— WJX+1/3+ ' V^^) 

For AdSd Schwarzschild, j = and joo = d— 1. (|13ip is undefined as both the denominator 
and the numerator vanish as j = 0. We should use the L'Hopital rule to obtain 

giujx—if)— ^g—iujx+if}— 



(132) 
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The vanishing of the denominator gives the locations of the QNM at cj. 



n 




6 = j3 + + The residue is also easily obtained as 



9± 



) 



sm 



TT(rf-l) 



res 



(133) 



Expanding f)12T[) as x — )• xo(r — > oo), we find the residue of the retarded correlator at QNM 
uj n given by 



which in the limit d — > 5 reduces to — 
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